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ABSTRACT 

We apply a simple model, tested on local ULIRGs, to disentangle the active galac- 
tic nucleus (AGN) and starburst contributions in submillimiter and 24 /im-selected 
ULIRGs observed with the Spitzer-IRS spectrometer. We quantitatively estimate the 
average AGN contribution to the stacked 6—8 iim rest- frame spectra of these sources in 
different luminosity and redshift ranges, and, under the assumption of similar infrared- 
to-bolonietric ratios as in local ULIRGs, the relative AGN/starburst contributions to 
the total infrared luminosity. Though the starburst component is always dominant in 
submillimeter-selected ULIRGs, we find a significant increase of the AGN contribution 
at redshift z > 2.3 with respect to lower z objects. Finally, we quantitatively confirm 
that the mid-infrared emission of 24 /im-selected ULIRGs is dominated by the AGN 
component, but the starburst component contributes significantly to the bolometric 
luminosity. 
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1 INTRODUCTION 

Infrared number counts and the cosmic infrared background 
provide strong evide nce that Ultralumi n ous Infrared Galax- 
Sande rs fc Mirabell 



Soifer et al. 1984' 



ies (ULIRGs; e .., 

11999 : iBlain et al.l |2002| : iLonsdale'eTal, .2006.). i.e. sources 
with a total infrared luminosity, I/8-ioooMm, greater than 
lO^^ Lfl, are a key component of the high redshift Universe 
e.g. iDole et all [200 ll : iLagache et al]|2005l : iLe Floc'h et~al] 



from a Spitzer GO-1 program to obtain low-resolution, mid- 
infrared spectra in the Spitzer First Look Survey (FLS). 
These 24 /tm-selected ULIRGs (hereafte r, 24 /^m-ULIRGs) 
have also redshifts in the range z = 1 — 3 (|Houck et al.ll2005l : 



I2OO5IV ULIRGs at high redshift have been discovered us- 
ing different selections, leading to distinct populations. In 
particular, observations performed with the Submillimeter 
Com mon User Bolometer Array (SCUBA; [Holland et al.l 
19991) on the Jarn es Clerk Maxwell Telescope (JCMT; e.g. 



Smail et all 1 19971 ) discovered a new population of submil- 
limeter galaxies (SM Gs), distributed at z = 1 — 3 (median 
redshift is 2.3) (e.g., IChapman et akllioosl ). and with esti- 
mated total infrared luminosity in the ULIRG range. We 
refer to these galaxies as submm-ULIRGs. 

Another population of high-z infrared luminous galax- 
ies has been discovered throu gh 24 /xm sel e ction. For exam- 
ple, iHouck et ai] (|2005l ) and iBrand et all (|2007l ) have sur- 
veyed the Bootes field of th e NOAO Deep Wide-Field Survey 
(NDWFS; Ijannuzi fc DevllT999.) wi th the Multiband Imag- 
ing Photo meter for Spitzer (MIPS; iRieke et al.|[2004l ) at 24 
fim. Also. FVan et al.l ( 20051 ) have reported the initial results 
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lYan et al.|[2005l : lBrand et al.lboO'iD . 

The properties of these high-z galaxies suggest that they 
are likely associated with an early phase in the formation of 
massive galaxies, supermassive black holes, and active galac- 
tic nuclei (AGNs). Therefore, making clear the origin of their 
energy sources for different redshift ranges is very important 
to understand the star formation history and the AGN ac- 
tivity in the high-z Universe. In this context, a key point 
to be investigated is the relative AGN/starburst contribu- 
tion to the bolometric luminosity of these galaxies. Here we 
present a quantitative estimate of the AGN content in sub- 
millimiter and 24 /xm-ULIRGs, based on a deconvolution 
method recently dev eloped by our group for local ULIRGs 
(jNardini et all (20081 . hereafte r N08) by means of Spitzer- 
Infrared Spectrograph (IRS; iHouck et alll2004l ) 5 - 8 fim 
spectroscopy. Throughout this paper we use a standard cos- 
mology with Ho ~ 73kms~^Mpc~^, SIm ~ 0.3 and ~ 0.7. 



2 SAMPLE DATA AND DATA REDUCTION 

Our sample consists of all the submm- and 24 /im- 
selected sources with known redshift, and with an avail- 
able Spitzer-IRS observation. The submm sample was ob- 
tained by combining the SCUBA detected sources from 
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Table 1. Submm- and 24 /xm-ULIR Gs sample. Column 1: Ob ject position. Col. 2, 3 , an d 4: Redshift and in frare d luminosity, in 
units of 1 0^3 Lp) (t jPope et al.ll200^') and re f erenc es th erein; [1] IChaprnan et al.l ||2005^: [21 IValiante et al.l l|2007h : [3] ISwinbank et al.l 
||2004| '): [4l lGreve et al.l | |2005^ : fSl lHouck et al.l | |2005|) : f6l lYan et al.l | |2005^ : f7l lBrand et alt ||2007| ') . Col. 5: Programs of the Spitzer-IRS 
observations. 
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Figure 1.6 — 8 fim rest-frame composite spectra of submm-ULIRGs and 24 ^m-ULIRGs, with their best fit models. The flux density, 
Fi, are normalized at 6.5 fim. Left: Composite spectra for all submm-ULIRGs (solid blue line) and 24 ^im-ULIRGs (dashed red line). 
Middle and right: Composite spectrum of the submm-ULIRGs and 24 /im-ULIRGs, respectively, with the best fit model {dotted line) 
and the AGN {dashed dotted line) and starburst (dashed line) components, and ctbol the intrinsic AGN contribution to the 6 fira 
flux density (in percent) and the AGN contribution to the total infrared luminosity (in percent), respectively, rg is the optical depth to 
the AGN at 6 ij,m. 



IChapman et al. (|2005|) (subm m bright, pre-selected in ra- 
dio) and lValiante et al.l l|2007l ) (without radio pre-selection) . 
These two subsamples are different in several respects (lens- 
ing properties, submm flux range, redshift distribution). 
However, their selection criteria are not expected to intro- 
duce any bias in their AGN content (we discuss possible ef- 
fects in §4). The 24 fim sources are also selected in different 
ways through their 24 ^m flux and different optical/infrared 



colors (e.g. 24/8 ^im , 24 /im/R). It is well known that these 
criteria are strongly biased toward sources whose mid-IR 
emission is dominated by the AGN component. Redshift, 
bolometric luminosities as well as other relevant data for 
all the sources are presented in Table 1. The bolometric lu- 
minosities listed in Table 1 are obtained by fltting a spec- 
tral energy distribution to the available data (see single ref- 
erences for details). The uncertainties on the single values 
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Figure 2. Upper panel: Absorption-corrected 6 /im/bolometric 
ratio versus the AGN fraction, ag, for submm-ULIRGs {blue 
filled circle) and 24 /im-ULIRGs sources {red open triangle). The 
continuous line refers to the best fit relation obtained for local 
ULIRGs (N08). Typical fit errors of x- and y- coordinates are a 
factor of 10% and 2 — 3, respectively. Lower panel: same for the 
stacked spectra of the subgroups discussed in the text. 



can therefore be high (about 50% for submm sources, and 
a factor of three for 24 fim sources, since the mid-IR to 
bolometric correction is strongly dependent on the AGN/SB 
relative contribution). Data reduction has been performed 
with the SPICE 1.4.1 package, using the pipeline version 
S15.3, provided by the Spitzer Science Center. The coad- 
ded images have been background-subtracted by means of 
the two observations in the nodding cycle. The calibrated 
one-dimensional spectra for the positive and the negative 
beams were extracted using the "optimal" extraction mode 
in SPICE, and the two one-dimensional spectra averaged in 
order to obtain the final spectrum. Since we are interested 
in the average properties of the two populations, and the 
spectral quality is low, we stacked the spectra of the single 
sources, after shifting and re-sampling each spectrum to the 
rest-frame wavelengths and bin sizes. We also had a quick 
look at each single spectrum, in order to check that the 
stacked data are not dominated by individual objects with 
high S/N. 



3 AGN AND STARBURST DECOMPOSITION 
MODEL 

In order to separate the AGN and the starburst components 
for submm-ULIRGs and 24 /^m-ULIRGs, we use the spectral 
decomposition model for local ULIRGs of N08. The model 
is based on the following main points: 1) in ULIRGs, the 
ratio between the 3 — 8 fim emission and the bolometric 
flux is 30 — 100 times higher for AGNs than for starbursts. 
The contrast rapidly decreases with increasing wavelength. 
2) The starburst spectrum in this wavelength range shows 
little variation from source to source, and a template ob- 
tained by averaging the highest signal-to-noise pure star- 
burst spectra fits well all the known starburst-dominated 
ULIRGs. 3) The AGN componen t can be repr esented with 
a power law f^'^^ cx A°-^ (Nctz er et al.ll2007l 'l. and an ex- 
ponential attenuation exp(— r(A)), where t he optical de pth 
follows the conventional law r(A) oc A~^'^^ jDrainelll989l ). In 
this assumption, the total flux density is written as follows. 



robs/ix rSB I J. AGN ~ 



r(A) 



(1) 



where rjsB and 77AGN are the flux density amplitudes of the 
starburst and AGN templates /^^ and /^*^^ normalized at 
6 /im, respectively. We can estimate the only two free pa- 
rameters (the ratio of r;sB to t^agn and re, which is the 6 
^m optical depth to the AGN) by fitting the spectrum. From 
these, we can get the intrinsic AGN contribution to the 6 
^m flux density, ae = »?agn/(?7agn + Vsb)- In addition, we 
can also estimate the AGN contribution to the total infrared 
luminosity (this roughly corresponds to the bolometric lu- 
minosity, Lboi, for ULIRGs), aboi = t?agn/(?7agn + Kriss), 
where K = R^'^^/R^^, and R^'^^ and R^^ are the ratios of 
the intrinsic flux at 6 to the total infrared flux of AGN 
and starbursts, respectively; in local ULIRGs, log R^'^^ = 



logi?^ 



1.911^:^2: yielding 7^-35 (Nardini 



et al. in prep.). The model has been successfully applied to 
a sample of local ULIRGs, obtaining two important results: 
1) despite the complexity of the spectra, all the sources were 
successfully fitted, showing that the relative AGN/starburst 
contribution and the extinction of the AGN component are 
responsible for most of the observed variety; 2) the expected 
AGN/starburst contributions to the bolometric luminosity 
reproduce closely the total observed luminosity. This is an 
important test for our model: the bolometric contributions 
are obtained from the 6 /xm spectral decomposition orfly 
(and from the average bolometric ratios, which have a fixed 
value for all objects), therefore the comparison between the 
predicted and observed total luminosities is an independent 
check of the results. 

Here we apply the same model to high-redshift sources. 
In doing so, we assume that a) the intrinsic SED, and b) the 
bolometric ratios, are the same at low and high redshift. We 
discuss the implications and the limits of these assumptions 
in U. 



4 RESULTS: AGN/SB DECONVOLUTION 
AND REDSHIFT DEPENDENCE OF AGN 
ACTIVITY 

The stacked spectra for the two samples of submm- and 
24 ^m-selected galaxies are shown in Fig. 1, together with 



4 Y. Watabe et al. 



our deconvolution in the AGN and starburst components. 
The main results are the following. 

• With only two free parameters, the deconvolution model 
reproduces the observed average spectral emission of high 
redshift sources in the 6-8 /xm range, similarly to what hap- 
pens for the low z ULIRGs analized in N08. We note that 
the simple fit we applied is strongly based on the assumption 
of little dispersion of the AGN and the starburst emission 
with respect to our templates. Therefore, it is not possible 
to extend the method to longer wavelengths (considering the 
redshift of our sources, we have a good spectral coverage up 
to 12 fim) due to the presence of extra features (mainly 
the silicate absorption) which are well known to have quite 
different strengths from source to source. 

• We quantitatively estimate the intrinsic AGN contribu- 
tions in submm-ULIRGs to the 6 fim and bolometric lumi- 
nosities: ae = 72.217'^%, Qboi = 6.9l^ g, where the errors 
refer only to the formal statistical uncertainties in the fits. 

• For the 24 /um-ULIRGs we find ae = 94.4l?;^% and 
Qboi = 32.5l3;g % . It has already b een suggested in pre- 
vious studies (e.g.. [Saiina et al]|2007l ). and confirmed from 
the AGN and starburst 5 — 8 lira to bolometric ratios (N08), 
that 24 ^m-ULIRGs are dominated by the AGN compo- 
nent. The new result presented here is that the starburst 
component is nevertheless present, and dominates the total 
bolometric luminosity. 

The above results rely on a series of assumptions which 
need to be tested, in order to assess the reliability of our 
method. In particular, as mentioned in §3, our choice of the 

5 — 8 /im rest-frame spectral interval is motivated by (1) the 
constancy of the AGN and starburst spectra, and (2) the 
high AGN/SB flux ratio in this band. These two properties 
are well measured for low-redshift sources, and the strongest 
assumption of our work is that they are both valid at higher 
redshifts and luminosities. Moreover, the possible intrinsic 
scatter in AGN and SB templates introduces additional un- 
certainties on top of the statistical errors in ag and Ofboi. A 
further uncertainty has to be added to the estimate of Qboi 
due to the possible intrinsic scatter in the AGN and the SB 
bolometric ratios. 

Following N08, it is possible to test the method and give 
a realistic estimate of the uncertainties by comparing the 
measured values of ae with the intrinsic 6 /xm/bolometric 
ratio R = 1/(6 /im)/Lboi. This new quantity has two inter- 
esting properties: 1) it depends only mildly from our spec- 
tral modeling (only through the AGN absorption rg) and 
never uses our assumptions on the bolometric ratios; 2) it 
is strongly sensitive to the AGN content, being (at low red- 
shift) ~ 35 times higher in AGNs than in starbursts. 

We plot R versus aa in Fig. 2. The most relevant re- 
sults are: (1) the two quantities are correlated as expected, 
i.e. the sources with the highest ae also have the highest 

6 /.im/bolometric ratio; (2) the continuous line in Fig. 2 
represents the best-fit correlation obtained by N08 for lo- 
cal ULIRGs, and the overall good match between the curve 
and our points suggests that there are no strong differences 
in bolometric ratios between low- and high- redshift sources. 
The dispersion of the single points from the relation is of the 
order of a factor ~2, while the values for the composite spec- 
tra have a dispersion of ~ 10-20%. We consider these as the 
best estimates of the errors involved in our analysis, and a 
confirmation of the validity of our approach. 
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Table 2. Intrinsic AGN contribution to the 6 /xm flux density, 
061 and to the total infrared luminosity, Obol in each group of 
submm and 24 fim sources discussed in the text. See Fig. 4 for 
the definition of Reg. 1-4. 




Figure 3. Same as the left panel of Fig. \T\ for submm-ULIRGs 
at z<2.3 {dotted lines) at z>2.3 (dashed lines.) The flux density, 
are normalized at 6.5 fim. 



In order to further investigate the issue of the AGN/SB 
content of high redshift galaxies, we analyzed the composite 
spectra obtained from different redshift and infrared lumi- 
nosity intervals. 

In Fig. 3, we show the composite spectra for the submm 
sources at z < 2.3 and z > 2.3, res pectively (z = 2.3 , is th e 
submm-ULIRG median redshift of IChapman erahl l|2005h . 
and roughly divides our sample in two subsamples of sim- 
ilar size). It is apparent from Fig. 3 that the polycyclic 
aromatic hydrocarbon (PAH) emission features of lower-2 
(z<2.3) submm-ULIRGs are stronger than those of higher-2 
{z>2.3) galaxies. The AGN contribution is not needed in the 
fit of the low-z sample, while (a6[%], Qboi[%]) = (87.2+3;^, 
I6.3I3-2) for the high-z sample. 

This correlation needs to be checked for several selec- 
tion/physical effects which could affect its significance and 
interpretation. In particular, we considered possible effects 
due to (1) selection, (2) K-correction , and (3) luminosity: 
1) Selection eff ects. The submm galax ies in our sample are 
obtain ed from I Chapman et al.l (|2005l ) and IValiante et al.l 
(|2007l ). The different selection criteria, based on radio pre- 
selection, lensing properties and submm flux are not ex- 
pected to be relevant for our results. However, we note that 
almost all the lValiante et ah ||2007|) objects have z>2.3 (Ta- 
ble 1), while those of lChapman et al.l (|2005h span the whole 
redshift range. In order to check whether the redshift effects 
in Fig. 3 are related to the properties of the two different 



AGN /starhurst 



content in high-z ULIRGs 5 




6 6.5 7 7.5 8 6 6.5 7 7.5 8 6 6.5 7 7.5 8 



Figure 5. Same as the middle panel of Fig. [T] for different subgroups of submm-ULIRGs: low-2, low luminosity sources (left), high-2;, 
low luminosity sources (middle), and liigh-z, high luminosity sources (right). 
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Figure 4. Luminosity-redshift plot for our samples of submm- 
ULIRGs {blue filled circle) and 24 /xm-ULIRGs {red open tri- 
angle). We define four regions delimited by the boundaries at 
2 = 2.3, LiR = 0.4 X 10" Lq and Ljr = IO^Lq. 



sampl es, we repeated our analysis for the IChapman et al.l 
(|2005l ) sample only, obtaining results perfectly compati- 
ble with those in Fig. 3. In o ther words, no signi ficant 
difference is found be tween the IValiante et al.l l|2007l ) and 
IChapman et al.l |200^ objects at z>2.3. 

2) K-correction effects. The SCUBA 850 /im selection does 
not seem to introduce any particular bias in favor of one 
of the two classes of sources; moving from redshift z — 1.8 
to z = 2.7 (the median redshifts for the two subsamples of 
submm-ULIRGs), the observed central wavelength changes 
from A = 300 /^m to A = 230 fim. Indeed, we checked that 
there is no difference between the AGN and starburst bolo- 
metric ratios at these two wavelengths by using the SEDs of 
the ULIRGs classified as AGNs and starburst, respectively 
l|Yanget al. | [2007 V 

3) Luminosity effects. Since the AGN fraction in lo- 
cal ULIRGs increas e s with increasin g infrared luminosity 
l|Veilleux et al.lll997l :l Kim et al.lll99d ). our results may im- 
ply that the real physical trend is with luminosity, which 
is obviously partially degenerate with redshift. In order 
to check for this possibility, we analyze separately the lu- 
minosity and redshift dependence. In Fig. 4 we show the 
luminosity-redshift plot of our sample. Analyzing the popu- 
lation of the four highlighted regions, we note that we can 



disentangle the luminosity and redshift effects by building 
composite spectra for submm-ULIRGs in regions 1 (L > 
10" Lq,z > 2.3), 3 (L = 0.4 - 1 X lO" Lq,z < 2.3) and 
4 (L = 0.4 - 1 X 10" Lq,z > 2.3). The results of the fits 
of these templates are presented in Table 2. Objects in re- 
gions 1 and 4 are in the same redshift interval but have quite 
different luminosities, while sources in regions 3 and 4 have 
similar luminosities but different redshifts. 

The spectra and the results are shown in Fig. 5. For 
submm-ULIRGs, the redshift dependence is confirmed by 
the comparison between the average spectra from the same 
luminosity range, while no strong luminosity dependence is 
found by comparing the composite spectra ai z > 2.3. The 
increase of the AGN fraction with the redshift, and its con- 
stancy with the luminosity are also indicated by the average 
6 /im/bolometric ratios, plotted in Fig. 2 as a function of 
the average AGN contributions. 

We have repeated the analysis the L — z plane also for 
the 24 fim sources, building compositte spectra for Reg. 1, 2 
and 3. The fit results presented in Table 2 indicate that no 
luminosity or redshift effect can be significantly detected. 



5 CONCLUSIONS 

We quantitatively estimated the relative AGN and star- 
burst contributions in two samples of submilliter- and 24 
/im-selected galaxies in the redshift range z — 1 — 3 hy ap- 
plying a deconvolution model based on AGN and starburst 
templates, already tested on low-redshift sources to the the 
stacked 6 — 8 /im spectra. We showed that the method gives 
consistent results on high redshift sources, suggesting that 
both the average 6 — 8 /im spectra and the 6 /im/bolometric 
ratios are similar in local and high redshift ULIRGs. Overall, 
we find that the starburst component is dominant in submm- 
ULIRGs both at 6 /im and bolometrically. We find evidence 
of a higher AGN contribution in higher redshift {z > 2.3) 
submm-ULIRGs with respect to lower redshift (1 < 2 < 2.3) 
objects. The AGN component always dominates the mid- 
infrared emission in 24 /xm-ULIRGs, as already known from 
previous work. However, the starburst contribution to the 
total luminosity is of the same order, or slightly higher, than 
that of the AGN even in the 24 /tm-ULIRGs. 
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